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ABSTRACT 


The  effects  of  roll  on  the  free-flight  yawing  motion  of  statically  stable 
test  models,  having  linear  force  and  moment  characteristics,  are  dis- 
cussed  in  more  detail  than  in  previous  publications.  The  correspondence 
between  the  damping,  frequency,  and  amplitude  of  the  nutations!  and  pre¬ 
cessions!  vectors  which  define  the  motion  is  listed.  Limitations  in  the 
use  of  simplified  data  reduction  procedures  when  the  model  motion  is 
monitored  in  only  one  plane  are  pointed  out,  and  cases  of  pseudo-nonlinear 
characteristics  being  exhibited  by  the  orthogonal  components  of  the  yawing 
motion  of  a  rolling  motion  are  demonstrated. 
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NOMENCLATURE 


CD 


Drag  coefficient 


C 

mPa 

cmq  +  Cm^ 


Static -st ability  derivative 
Magnus -moment  derivative 

^  C  d  G 

Damping- in-pitch  derivatives,  - a —  +  - — !=— 

(i/v)  da  (i/v) 


CN„ 

H,  M,  P,  T 

Iv 


Kl,K2 

ki,k2 

ka 

kt 

S. 

m 

P 

S 

V 

x 

a, 


6 

01>  &2 
Pi'  P2 
5 

P 

/  f 

<i>l>  $2 


Normal-force  derivative 
Constants  in  Eq.  (1) 

Model  moment  of  inertia  (relative  to  the  longitudinal 
axis) 

Model  moment  of  inertia  (relative  to  a  transverse 
axis) 

Constants  in  Eq.  (2) 

Absolute  values  of  Kx  and  K2,  respectively 
Radius  of  gyration  relative  to  the  longitudinal  axis /i 
Radius  of  gyration  relative  to  a  transverse  axis/i 
Reference  length 
Model  mass 

Model  roll  rate  (with  respect  to  distance  traveled) 
Reference  area 
Model  velocity 
Distance  traveled 

Components  of  the  complex  yaw  angle 

V/3J  +  a1 

Constants  in  Eq.  (2) 

Damping  rates  of  the  motion  vectors  in  Eq.  (2) 

Complex  yaw  angle 

Mass  density  of  the  range  air 

Rates  of  rotation  of  the  motion  vectors  in  Eq.  (2) 
Precession  angle,  see  Fig.  1 
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SUPERSCRIPTS 

/ 

t, 

SUBSCRIPTS 

o 

n 

P 


First  derivative  with  respect  to  distance 
Second  derivative  with  respect  to  distance 


Corresponding  to  the  start  of  a  motion  pattern 
Nutational  vector 
Precessional  vector 
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SECTION  I 
INTRODUCTION 


In  free -flight  testing  in  either  an  aeroballistic  range  or  in  a  wind 
tunnel,  certain  stability  derivatives  can  be  evaluated  from  the  measured 
yawing  motion  of  the  model.  For  a  nonrolling,  axisymmetric  model,  the 
damping  parameter  and  the  frequency  of  each  of  the  two  orthogonal  com¬ 
ponents  of  the  yawing  motion  (o  and  ,3)  are  constant;  further,  both  com¬ 
ponents  have  identical  damping  parameters  and  identical  frequencies. 

As  the  roll  velocity  of  the  model  increases  from  zero,  the  complexity 
of  the  motion  of  each  orthogonal  component  increases  and  consists  of 
two  subcomponents  having  different  damping  parameters  and  different 
frequencies.  It  is  apparent  from  the  above  that  the  analysis  of  the  yaw¬ 
ing  motion  of  a  rolling  model  can  be  appreciably  more  complicated. 

Considerable  emphasis  is  currently  being  directed  to  the  free-flight 
testing  (both  aeroballistic  ranges  and  wind  tunnels)  of  aerodynamic  con¬ 
figurations  that  are  statically  stable  but  subject  to  experience  small  roll 
rates.  The  small  roll  rates  arise  from  disturbances  during  the  launching 
of  the  free-flight  models,  from  small  asymmetries  in  nominally  symmetric 
models,  or  from  built-in  asymmetries;  this  is  in  contrast  to  the  large  roll 
rates  of  spin-stabilized  projectiles.  It  follows  that  a  basic  concern  in  free- 
flight  testing  (particularly  in  wind  tunnels  when  motion  is  viewed  in  one 
plane)  is  related  to  how  large  a  roll  velocity  can  exist  and  it  still  be 
assumed  zero  in  reducing  stability  data  from  the  measured  model  motion. 

The  purpose  of  this  report  is  to  present  a  discussion  of  the  effects  of 
roll  on  the  yawing  motion  of  free-flight  models  experiencing  small  roll 
rates.  In  particular,  the  significance  of  roll  effects  in  data  reduction 
procedures  is  considered. 


SECTION  II 

ANALYSIS  OF  FREE-FLIGHT  RANGE  DATA 


The  basic  roll  effects  discussed  in  this  report  correspond  to  axisym¬ 
metric  configurations.  In  discussing  the  effects  of  model  roll  on  the  yaw¬ 
ing  motion  of  a  test  model,  it  is  useful  to  examine  the  conventional  data 
analysis  procedure  as  used  in  aeroballistic  range  testing.  In  range  test¬ 
ing,  roll  rates  from  zero  to  large  values  consistent  with  spin- stabilized 
projectiles  are  experienced;  hence,  data  reduction  procedures  permitting 
a  nonzero  roll  velocity  are  necessarily  required. 
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In  range  work  the  differential  equation  normally  used  in  describing 
the  general  rolling-yawing  motion  of  an  axisymmetric  model  is 

f  +  (H  -  iP)  i  -  (M  +  iPT)  £  =  0  (!) 

where 

H  =  (pS/2m)  TCLa  -  CD  -  kt“3  ( Cmq  +  C^J 
M  =  (/>S/2m£)  kt  2  CfrlG 
T  =  (pS/2m)  £Cl0  +  ka  1  Cmp  J 
P  =  ( I  x  /I  y )  P 

Derivations  of  this  equation  (relative  to  a  nonrolling  axis  system)  have 
been  presented  previously,  for  example.  Refs,  1  and  2.  The  solution  of 
Eq.  (1)  for  a  model  having  linear  variations  of  force  and  moment  with 
yaw  angle,  and  a  constant  p  can  be  written 


€  = 


(2) 


where 

di  "  ^  + 


The  0's  are  roots  of  the  auxiliary  equation  corresponding  to  the  differ¬ 
ential  equation  designated  Eq,  (1). 


It  should  be  noted  that  the  damping  of  the  motion  corresponding  to 
the  0j  root  is  defined  by  the  real  part  of  the  root  and  that  the  frequency  of 
the  motion  is  defined  by  the  imaginary  part  of  the  root.  Hence,  for  oscil¬ 
latory  motion  to  exist  the  root  must  have  a  nonzero  imaginary  part. 

The  roots  expressed  as  functions  of  the  coefficients  of  Eq.  (1)  are 


ii,  +  i  tf>t  =  1/2  [-H  +  iP  +  y'4M  4-  H2  -  P2  +  i  2P  (2T  -  H)  ]  (3) 

and 

fx2  +  i  <f>'2  =  1/2  [-  H  +  iP  -  V'lM  <  H1-  P*  +“  i  2P  (2T  -  H)]  (4) 

Values  for  the  K's  and  0's  of  Eq.  (2)  are  obtained  by  fitting  Eq.  (2) 
to  the  measured  variations  of  the  ct  and  /3  components  of  the  angular 
motion  of  the  model  with  distance  traveled.  This  equation  is  nonlinear 
in  terms  of  the  0's,  and  the  use  of  an  iterative,  differential  correction- 
type  curve -fitting  procedure  is  necessary.  The  desired  H,  M,  P,  and 
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T  parameters  can  then  be  determined  as  functions  of  the  evaluated  6fs 
with  use  of  the  expressions  for  the  sum  and  the  product  of  (p^  +  and 
Cm 2  +  i$2)  and  can  be  written 

H  -  -[p.  +  (5) 

P  =  <^>1  +  4*2  (6) 

T - !  p  (7) 


M  =  4>\4>\  -  v-iHi 

*  * 

“  ^1  5>1 


(8) 


The  assumption  in  the  approximate  relationship  for  M,  that  the  produc 
of  the  damping  parameters  is  small  relative  to  the  produce  of  the  frequen¬ 
cies,  is  quite  reasonable  for  a  typical  aerodynamic  configuration. 

The  above  procedure,  permitting  large  roll  velocities,  is  very  ade¬ 
quate  in  the  analysis  of  yawing  motion.  It  is  difficult  for  one  to  appreciate 
the  significance  of  the  roll  effects  that  the  procedure  accounts  for;  however 
considerable  insight  to  the  problem  of  roll  effects  can  be  obtained  by  exam¬ 
ining  the  related  equations.  If  P  is  set  equal  to  zero^  (implying  a  zero  roll 
rate)  in  Eqs.  (3)  and  (4),  the  corresponding  p 's  and  $'s  are 

H  +  i  4>'  -  1/2  [-  H  ±  \I  4\1  +  H2  ]  (9) 


Considering  that  oscillatory  motion  requires  a  complex  or  imaginary 
root  and  that  H  is  a  real  number,  the  term  (4M  +  H^)  must  be  less  than 
zero,  or  4M  <  -  H*.  From  the  definition  of  M  listed  previously, 

4M  =  4  (pS/2m£)  kf3  Cma 

-  and  it  follows  that  Cm  must  be  negative.  An  aerodynamic  configuration 
is  defined  as  being  statically  stable  when  it  has  a  negative  Cm  value. 
Hence  for  oscillatory  motion  (nonrolling  model)  it  is  apparent  from  Eq.  (9) 
that  the  model  must  be  statically  stable. 

The  two  terms  on  the  right  side  of  Eq.  (2)  defining  the  complex  yaw 
angle,  §,  correspond  to  rotating  vectors  in  the  {a,  plane  (see  Fig.  1). 
The  model  yawing  motion  defined  by  the  two  vectors  is  obviously  depend¬ 
ent  on  the  values  of  the  K's  and  0's.  Values  of  the  0's  are  functions 
primarily  of  the  aerodynamic  and  inertia  characteristics  of  the  model  and 
the  rolling  velocity;  the  K's  are  dependent  on  the  initial  disturbances  of 
the  model.  For  general  rolling-yawing  motion,  Kj,  K2,  6>i,  and  8 2  are 
general  complex  numbers;  however,  roll  effects  on  the  yawing  motion  of 
a  model  become  more  apparent  from  first  examining  the  more  restricted 
types  of  model  motion. 
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Fig.  1  Computed  Motion  Patterns  as  Functions  of  Initial  Disturbances  and  Model  Roll 

The  basic  types  of  model  motion  corresponding  to  assigned  restrictions 
on  the  K's  and  0's  of  Eq.  (2)  are  demonstrated  in  the  following  examples, 
in  which  zero  damping  is  assumed: 

a.  If  Ki  and  K2,  and  0 1  and  0g  are  complex  conjugates  then 
the  motion  is  planar  and  along  the  real  axis.  This  type 
corresponds  to  the  angular  motion  of  a  one-degree-of- 
freedom  dynamic  balance  system  as  used  in  a  wind 
tunnel.  A  basic  point  to  be  noted  here  is  that  for 

the  0's  to  be  complex  conjugates,  the  roll  velocity,  p, 
must  necessarily  be  zero,  which  is  apparent  from 
■  Eqs.  (3)  and  (4).  It  should  be  noted  that  the  restriction 
on  the  0's  dictates  that  the  angular  frequencies, 
are  equal  in  magnitude  and  have  opposite  signs. 

b.  Let  the  restriction  of  fii  and  02  being  complex  conjugates 
be  retained  {p  =  0).  If  K  \  and  K2  are  not  required  to  be 
complex  conjugates  but  Iki|  =  I K 2;  ,  then  general  planar 
motion  exists.  This  is  indicated  in  the  sketch  in  Fig.  la 
where  the  yawing  motion  is  along  a  line  displaced  from  the 
a  and  j3  axes. 

c.  If  the  restrictions  of  item  (b)  are  retained  except  that 

|  Ki  |  f  !  K2  |,  then  elliptic  motion  exists  as  indicated  in 
Fig  lb.  The  limiting  case  here  is  circular  motion  corre¬ 
sponding  to  one  of  the  K's  being  zero. 
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d.  If  Ki  and  K2  are  general  complex  numbers  as  in  item  (c) 
and  the  restriction  of  and  9 2  being  complex  conjugates 
is  removed,  the  resulting  motion  is  of  the  precessing 
elliptic  type  shown  in  Fig.  lc. 

Removing  the  restriction  of  Oj  and  9 2  being  complex  conjugates  corre¬ 
sponds  to  a  nonzero  roll  velocity  (see  Eqs.  (3)  and  (4)).  Hence,  the  pre¬ 
cessing  elliptic  motion  of  item  (d)  is  caused  by  model  roll,  and  the  motion 
patterns  of  items  (a,  b,  c)  are  characteristic  of  nonrolling  yawing  motion. 


The  relationship  between  the  roll  velocity  of  a  model  and  the  direction 
and  magnitude  of  the  precession  of  the  corresponding  elliptic  motion  is 
obtained  in  the  following  derivation.  In  Fig.  lc  the  angle  of  precession, 
i!j,  is  defined  as  the  angle  between  two  adjacent  peaks  as  indicated  in  the 
sketch.  The  peak  at  (a)  corresponds  to  a  point  in  flight  where  both  vectors 
are  aligned,  and  at  peak  (b)  the  vectors  are  aligned  again  after  each  vector 
has  rotated  approximately  360  deg.  As  previously  noted,  the  two  vectors 
rotate  in  opposite  directions  and  as  f-  | <6 2 1  ,  for  (p  f  o),  then  \p  must 
have  a  nonzero  value.  The  sign  of  ip  is  defined  consistent  with  rotating 
vector  notation;  hence,  if  the  motion  precesses  counterclockwise  ip  is 
positive.  It  is  apparent  from  the  above  that  the  motion  will  precess  in 
the  direction  of  the  rotation  of  the  larger  absolute  frequency.  From 
Eq.  (6)  it  follows  that  the  larger  absolute  frequency  will  have  the  same 
sign  as  the  roll  velocity  of  the  model.  Consistent  with  the  notation  of 
Fig.  lc, 

( 360°  +  $)  /  <j>[  =  (-360°  +  t>)  1  <b[ 

for 

1/2  <  |  4>\  l  4>\  I  <1 

which  corresponds  to  the  roll  rates  of  interest  in  this  report.  Hence, 

(<jh/&)  =  U  -  360°)  /  (ip  +  360°)  (10) 

With  use  of  Eqs.  (6)  and  (10)  an  expression  for  P  can  be  written  as 

P  =  [2 P  ;  {iif  +  360  )] 


Hence, 

p  =  (  ly  / 1  x  )  [  2  \jj  /  ( p  +  360 )  ]  (j)2 

(11) 

and 

(12) 

I P  =  360°  P  /  [  2  ( I  y  / 1  x )  -  p] 

For  these  expressions,  ip  is  in  degrees;  hence,  p  and  the  <£'s  are  in  deg/ft 
units. 
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Equations  (11)  and  (12)  show  explicitly  the  correspondence  between 
the  measured  precessing  elliptic  motion  of  a  test  model,  the  inertia  ratio 
(Iy/Ix),  and  the  roll  velocity  of  the  model. 

A  further  point  of  interest  related  to  an  (ft,  /3)  plot  is  roll  resonance 
which  can  be  of  concern  in  testing  an  asymmetric  configuration.  Roll 
resonance  occurs  when  p  =  <f> 2,  considering  that  the  condition  of  resonance 
is  slightly  dependent  on  the  damping  of  the  system.  The  precession  angle 
corresponding  to  resonance,  i£res,  can  be  obtained  from  Eq.  (12): 

^res  *  [  2(ly/]x)  _  1  ] 

In  the  above  examples  of  model  motion  defined  by  Eq.  (2),  it  should 
be  noted  that  when  |ki  J  t  |K2  J  ,  the  motion  sketches  are  for  the  case 
where  the  vector  having  the  larger  amplitude  corresponds  to  the  larger 
absolute  frequency.  The  correspondence  between  the  magnitudes  of  fre¬ 
quency,  damping,  and  amplitude  of  the  vectors  is  discussed  in  detail  in  a 
later  section  of  this  report. 


SECTION  111 

SIMPLIFIED  DATA  ANALYSIS 


As  previously  noted,  the  components  0  and  a  of  the  yawing  motion  of 
a  nonrolling  model  are  defined  by  damped  sinusoidal  curves  having  identi¬ 
cal  frequencies  and  damping  parameters.  It  follows  that  in  free-flight 
testing  in  wind  tunnels,  where  the  model  motion  is  monitored  in  only  one 
plane,  it  is  particularly  significant  to  be  able  to  assume  that  the  model  has 
nonrolling  motion.  However,  some  model  roll  arising  from  the  launching 
of  the  model  or  from  small  model  asymmetries  that  are  either  accidental 
or  designed  into  the  model  may  exist,  and  the  concern  in  this  type  of  test¬ 
ing  is  the  sensitivity  of  the  components  of  the  yawing  motion  to  model  roll. 
To  aid  in  examining  this  problem,  the  3  and  a  components  of  Eq.  (2)  can 
be  separated  and  written  as 

ft  =  (  a  cos  0,  x  -  b  sin  9^  x)  exp  +  (ccos<£,x  -  d  sin  <£,  x)  exp  (/i,  x)  (13) 

and 

a  ~  (bcoH^>[X  +  a  sin  ^>,  x)  exp  (fi[  x)  +  (  d  cos  ^  1  +  c  sin  <£3  x)  exp  x  )  (14) 

where  a,  b,  c,  and  d  are  real  constants. 

From  Eqs.  (13)  and  (14)  it  is  apparent  that  each  component,  for  the 
case  of  a  rolling  model,  is  the  sum  of  trigonometric  terms  containing 
different  frequencies;  hence,  the  resulting  curve  for  the  simplified  case 
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of  zero  damping  will  necessarily  have  an  apparently  varying  frequency 
and  a  varying  amplitude.  Examples  of  such  motion  are  shown  in  Figs.  2, 
3,  and  4.  The  motion  plots  can  be  deceiving  in  that  apparent  nonlineari¬ 
ties  can  be  indicated.  For  large  roll  rates,  Fig.  2,  these  apparent  non¬ 
linear  characteristics  in  a  component  of  the  motion  are  easily  detected; 
hence,  the  primary  concern  in  experiments  is  related  to  possible  effects 
on  the  yawing  motion  existing  at  small  roll  rates  where  the  apparent  non- 
linearities  are  not  obvious  from  the  motion  plots  when  the  model  motion 
is  monitored  in  only  one  plane.  The  model  motion  plots  in  both  Figs.  3 
and  4  correspond  to  zero  damping  and  a  precession  angle  of  5  deg  and  are 
representative  of  the  case  where  each  individual  component  of  the  motion 
appears  to  be  a  sinusoidal  trace  with  a  constant  damping  factor. 


-8  -6  -4  -2  0  2  4  6  8 

|3,  deg 

a.  (Alpha,  Beta)  Variation 

Fig.  2  Computed  Motion  for  a  Statically  Stable  Configuration,  il/  -  26.6  deg 
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As  can  be  observed  from  the  corresponding  (o,  |3)  plot,  in  Fig.  3, 
the  initial  disturbance  of  the  oscillatory  motion  {?0)  is  in  a  plane  oriented 
at  4  5  deg  relative  to  the  n  and  /3  axes,  and  in  Fig.  4,  ?0  i-s  i-n  the  plane 
corresponding  to  the  a  axis.  Precession  angles  well  above  5  deg  have 
been  measured  on  axisymmetric  models  launched  in  aeroballistic  range? 
using  unrifled  launch  tubes;  hence,  precession  angles  on  the  order  of 
5  deg  could  be  expected  in  wind  tunnel  tests,  particularly  with  asym¬ 
metric  test  models. 


-6  -4  -2  0  2  4  6 

P.  deg 

o.  (Alpha,  Beta)  Variation 

Fig.  3  Computed  Motion  for  a  Statically  Stable  Configuration,  &  =  5  deg, 
£o  Displaced  45  deg  from  the  a-Axis 
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The  variation  of  the  a  component  in  Fig.  4  is  more  representative  of 
the  case  confronted  in  wind  tunnel  testing;  this  follows  from  noting  that  the 
initial  angle  disturbance  in  wind  tunnel  tests  can  be  primarily  restricted  to 
the  plane  in  which  the  yawing  motion  is  being  monitored.  In  examining  the 
a  component  (Fig.  4),  presuming  roll  effects  are  negligible,  one  obtains  a 
decay  of  the  a  component  in  four  cycles  of  about  0.  4  deg,  although  the  com¬ 
puted  model  motion  corresponds  to  zero  aerodynamic  damping.  An  ampli¬ 
tude  decay  error  of  0.4  deg  corresponds  to  an  error  of  about  twenty  percent 
in  the  measured  damping  factor  of  a  system,  presuming  that  the  damping 
of  the  system  causes  an  amplitude  decay  in  the  four-cycle  interval  of  about 
0.3  of  the  initial  amplitude  of  7  deg,  which  is  quite  reasonable.  There  was 
no  measureable  difference  in  the  frequency  of  the  a  component  of  the  motion 
of  Fig.  4  from  the  exact  nonrolling  frequency. 


-6  -4  -2  0  2  4  6 

0,  deg 

a.  (Alpha,  Bela)  Variation 

Fig.  4  Computed  Motion  for  a  Statically  Stable  Configuration,  $  =  5  deg, 
£o  Aligned  with  the  c-Axis 
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The  above  analysis  indicates  that  the  wind  tunnel  data  reduction  pro¬ 
cedure  in  which  only  one.  component  of  the  yawing  motion  is  examined 
on  the  basis  of  presumed  negligible  roll  effects  in  some  cases  can  lead 
to  significant  errors  in  the  measured  damping;  hence,  particular  care 
should  be  exercised  in  this  type  of  an  analysis  to  ensure  that  roll  effects 
are,  in  fact,  negligible. 

It  should  be  noted  that  other  approximate  procedures  propused  for 
the  analysis  of  free-flight  data  and  dependent  on  measured  motion  decay 
(e.  g.  ,  see  Ref.  3)  must  be  used  cautiously;  this  will  become  more 
evident  in  a  later  section  of  this  report.  Further,  a  point  of  particular 
,  interest  concerns  the  curve  fitting  of  typical  motion  patterns  in  the  a-/ 3 
plane.  In  the  discussed  range  procedure,  the  measured  values  of  both 
a  and  /3  are  fitted  simultaneously,  and  hence,  various  types  of  motion 
from  planar  to  patterns  exhibiting  large  precession  can  be  fitted  ade¬ 
quately.  It  was  noted  in  Ref.  4  that  all  of  the  equation  constants  in  the 
expression  for  the  complex  yaw  angle,  are  present  in  each  component 
of  (a  or  3)-  Further,  it  was  stated,  in  view  of  this,  that  the  unknown 
equation  constants  could  be  determined  by  fitting  either  the  a  or  |3  values 
separately  or  both  simultaneously.  That  the  equation  constants  corre¬ 
sponding  to  general  motion  patterns  can  be  evaluated  adequately  by  fitting 
one  component  of  the  motion  is  believed  to  be  optimistic.  The  case  where 
the  motion  pattern  of  an  (o,  /3)  plot  can  be  adequately  fitted  using  only  a 
values,  for  example,  implies  that  the  motion  pattern  is  uniquely  defined 
by  only  the  a  values.  This  condition  can  exist  for  model  motion  for  which 
large  precession  angles  exist  (usually  associated  with  large  roll  rates); 
for  planar  motion,  it  is  apparent  that  the  motion  pattern  is  not  uniquely 
defined  by  only  one  of  the  components.  Hence,  there  will  be  a  transition 
region  between  planar  motion  (that  can  not  be  adequately  fitted)  and  motion 
having  large  precession,  and  in  this  region  difficulty  in  curve  fitting  can 
be  expected.  Unfortunately,  a  large  portion  of  current  stability  testing 
falls  in  this  transition  region. 


SECTION  IV 

CORRESPONDENCE  BETWEEN  <£’s,  ji’s,  K's.  and  P 


There  are  pronounced  variations  that  can  exist  in  the  basic  precess- 
ing  elliptic  motion  patterns  shown  previously.  From  data  reduction  con¬ 
siderations,  recognition  of  these  variations  can  be  particularly  significant 
(e.g.,  see  Ref.  5)  as  an  aid  in  confirming  whether  or  not  the  measured 
motion  is  basically  linear.  The  possible  variations  of  the  yawing  motion 
of  a  rolling  body  are  associated  with  the  correspondence  between  the  0's, 
n' s,  and  K's  in  Eq.  (2).  The  two  vectors  in  Eq.  (2)  have  been  identified 
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in  the  past  as  the  nutational  vector  (having  the  larger  absolute  frequency) 
and  the  precessional  vector.  However,  from  previous  publications  it  is 
difficult  to  identify  which  damping  or  which  amplitude  (as  to  the  smaller  or 
larger  value)  corresponds  to  the  nutational  or  precessional  vector.  The 
purpose  of  the  following  discussion  is  to  clarify  the  correspondence  between 
the  p' s,  o's,  and  K's  and  to  serve  as  an  aid  in  recognizing  various  possible 
variations  in  the  motion  patterns. 


The  correspondence  between  the  p's  and  $'s  can  be  better  defined  by 
evaluating  the  complex  radical  in  Eqs.  (3)  and  (4).  In  Ref.  1,  the  corre¬ 
sponding  radical  was  evaluated  approximately  by  using  a  binomial  expan¬ 
sion.  The  approximate  equations  obtained  for  the  p's  and  $'s  in  Ref.  1 
indicate  that  the  larger  damping  value  will  always  correspond  to  the  nuta¬ 
tional  vector.  The  use  of  approximate  solutions  for  similar  complex 
radicals  has  continued  in  recent  publications  (e.  g.  ,  Ref.  6).  However, 
the  evaluation  of  such  a  complex  radical  can  be  obtained  precisely  with 
use  of  DeMoivre's  theorem  and  can  be  written  as 


a  +  ib 


+ 


b  +  a 


(15) 


The  sign  of  the  imaginary  term  is  positive  if  (b)  (the  imaginary  part  of  the 
radical -Ja  +  ib)  is  positive,  and  is  negative  if  (b)  is  negative.  From 
Eqs.  (3)  and  (4),  a  =  4M  +  H2  -  p2  and  b  =  2P  (2T  -  H). 


In  evaluating  the  complex  radical  of  Eqs.  (3)  and  (4),  it  should  be 
noted  that  the  term  (2T  -  H),  in  general,  can  be  expected  to  be  negative; 
however,  it  is  quite  possible  for  the  term  to  be  positive  for  certain  con¬ 
figurations.  It  follows  then  that  the  sign  of  (b)  will  be  the  same  as  the 
sign  of  P  for  (2T  -  H)  >  0  and  opposite  the  sign  of  P  for  (2T  -  H)  <  0. 

With  use  of  Eq.  (15),  the  following  expressions  can  be  written  for  the  p's 
and  <£'s: 


ii,  =  1/2 

<f>\  =  1/2 

ii2  =  1/2 

<f>[  =  1/2 


-H  +  ^ 


V(4m  +  H*  -  P*)*  +  [ 2P  (  2T  -  H)]*  +  (4m  +  H*  -  P2) 

2 


P  -  (±> 
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Consistent  with  above  comments,  for  the  case  (2T  -  H)  <  0  the  positive 
sign  preceding  the  radical  in  the  expressions  for  the  <£'s  is  used  for  p  <  0 
and  the  negative  sign  is  used  for  p  >  0.  The  above  expressions  indicate 
that  the  larger  damping  value  corresponds  to  the  nutational  frequency  when 
(2T  -  H)  <  0,  however,  for  configurations  where  {2T  -  H)  >  0,  the  larger 
damping  value  corresponds  to  the  precessional  frequency. 

Although  the  (±)  sign  on  the  right  side  of  Eq.  (15)  precedes  the  imagi¬ 
nary  term,  it  should  be  pointed  out  that  Eq.  (15)  would  hold  equally  as  well 
if  the  (±)  sign  preceded  the  real  term  and  was  used  consistent  with  the 
restriction  on  the  sign  of  (b).  Such  a  sign  change  would  not  affect  the  above 
comments  on  the  correspondence  of  the  ju's  and  $'s. 

The  relationship  between  the  frequency  of  a  vector  and  its  amplitude 
is  dependent  on  the  initial  disturbances  of  the  model  and  can  be  examined 
with  use  of  the  following  sketches  in  which  the  yawing  motion  has  the  same 
initial  displacement,  £0: 


Pattern  When  Nutational  Vector  Pattern  When  Precessional  Vector 

Has  the  Larger  Amplitude  Has  the  Larger  Amplitude 

Sketch  a  Sketch  b 


Note  that  the  ?  vector  in  Sketch  (a)  is  rotating  counterclockwise  and 
is  rotating  clockwise  in  Sketch  (b).  It  is  apparent  that  the  f  vector  will 
rotate  (for  small  p  values)  in  a  direction  consistent  with  the  rotation 
direction  of  the  motion  vector  having  the  larger  amplitude,  (Kn  or  Kp). 

As  previously  discussed,  the  elliptic  motion  will  precess  in  the  direction 
of  the  rotation  of  the  nutational  vector  (direction  of  the  model  roll);  hence, 
the  nutational  vector  will  have  the  larger  amplitude  only  if  the  £,  vector 
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rotation  is  in  the  same  direction  as  that  of  the  nutational  vector.  The 
direction  of  rotation  of  the  I  vector  and  the  ratio  of  the  amplitudes  of  the 
two  vectors  corresponding  to  the  elliptic  motion  is  defined  by  the  initial 
velocity  disturbance  of  the  ?  vector,  (|0). 

Some  summary  comments  concerning  the  nutational  and  precessional 
vectors  which  follow  from  the  above  discussion  of  precessing  elliptic 
motion  are; 

a.  The  angular  direction  of  precessing  elliptic  motion  will  be 
the  same  as  the  rotation  direction  of  the  nutational  vector; 
the  nutational  vector  will  always  have  the  same  rotation 
direction  as  that  of  the  model  roll. 

b.  The  nutational  vector  will  have  the  larger  damping  value 
if  (2T  -  H)  <  0,  and  conversely  if  (2T  -  H)  >  0. 

c.  The  nutational  vector  will  have  the  larger  amplitude  only 
if  the  rotation  of  the  f  vector  has  the  same  direction  as 
the  rotation  direction  of  the  nutational  vector.  The  direc¬ 
tion  of  the  rotation  of  the  f  vector  and  the  ratio  of  the 
amplitudes  of  the  two  motion  vectors  are^defined  by  the 
initial  velocity  disturbance  of  the  model,  Co- 

d.  In  consideration  of  item  (c),  it  follows  that  the  decay  of 
the  amplitude  peaks  on  an  (a,  £)  plot  is  dependent  on  the 
initial  disturbances  of  the  model. 

Figure  5  furnishes  an  indication  of  the  variations  in  the  damping 
parameters  and  the  frequencies  of  the  motion  vectors  as  functions  of  the 
model  roll  (precession  angle)  for  assigned  aerodynamic  and  inertia 
parameters  representative  of  a  5 -deg  semiangle  cone  model. 

From  the  correspondence  between  the  /Fs  and  0's  and  between  the 
$'s  and  the  amplitudes  of  the  vectors,  it  follows  that  a  model  can  be 
disturbed  such  that  the  larger  amplitude  vector  can  have  either  the 
larger  or  smaller  damping  value;  this  is  significant  in  that  it  implies 
that  the  variation  in  the  amplitude  of  the  yawing  motion  of  a  rolling 
model  is  a  function  of  the  initial  disturbances  of  the  model.  In  Ref.  3 
an  approximate  data  reduction  procedure  is  proposed  in  which  it  is 
assumed  that  a  curve  of  the  logarithm  of  the  absolute  value  of  C,  (6), 
corresponding  to  the  peaks  of  the  (a,  fi)  plot  as  a  function  of  distance 
traveled  is  linear  and  that  the  slope  of  this  curve  is  proportional  to  the 
damping  of  the  configuration.  Values  of  the  damping  factor  (jug)  evalu¬ 
ated  consistent  with  the  procedure  of  Ref.  3  have  been  obtained  from 
computed  model  motion  patterns  corresponding  to  a  5-deg  semiangle 
cone  configuration. 
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The  values  obtained  for  two  different  initial  disturbance  conditions 
are  shown  in  Fig.  6  and  indicate  that  the  damping  evaluated  as  a  function 
of  the  measured  amplitude  is  dependent  on  the  initial  disturbances  of  a 
rolling  model.  Deviations  in  the  measured  damping  can  be  appreciable 
at  the  larger  roll  velocities. 


Fig.  6  Damping  Evaluated  from  the  Amplitude  Variation  of  the  Computed  Motion  for  a 


Simulated  5-deg,  Semiangle  Cone  Model,  /i.j  > 


Note  from  the  correspondence  between  ^'s,  $'s,  andK's  that  the  £ 
vector  rotates  in  a  direction  consistent  with  the  rotation  of  the  motion 
vector  having  the  larger  amplitude  (for  small  roll  rates);  it  follows  when 
the  larger  amplitude  vector  initially  has  the  larger  damping,  the  motion 
(?  vector)  can  have  a  reversal.  This  type  of  (a,  j3>  plot  for  computed 
motion  is  shown  in  Fig.  7  where  the  precessional  vector,  kp,  initially 
has  the  larger  damping  and  the  slightly  larger  amplitude.  Later  in  the 
flight  the  precessional  vector  has  damped  such  that  it  has  the  smaller 
amplitude  and  a  motion  reversal  occurs. 
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The  motion  patterns  in  Fig.  7  and  Fig.  2  are  of  particular  significance 
in  conjunction  with  comments  made  in  Ref.  5.  In  Ref.  5  motion  patterns 
exhibiting  converging-diverging  oscillations  and  motion  reversals  were 
attributed  to  probable  nonlinearities  in  the  model  motion.  Although  the 
test  configurations  of  Ref.  5  may  have  had  nonlinear  characteristics,  it 
should  be  pointed  out,  in  consideration  of  the  above  discussion,  that  meas¬ 
ured  converging-diverging  oscillations  and  motion  reversals  are  not  neces¬ 
sarily  indicative  of  nonlinear  model  motion. 


0,  deg 

a.  (Alpha,  Beta)  Variation 

Fig.  7  Computed  Motion  tor  a  Simulated  Linear  Configuration  to  Demonstrate  Motion  Reversal 
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SECTION  V 

CONCLUDING  REMARKS 


The  effects  of  roll  on  the  free -flight  yawing  motion  of  statically  stable 
test  models  having  linear  force  and  moment  characteristics  are  discussed 
in  more  detail  than  in  previous  publications. 

The  zero-damping  motion  patterns  defined  by  the  basic  linear  equa¬ 
tion  of  motion  for  an  axisymmetrie  model  are  described  with  emphasis 
on  the  effects  of  model  roll  and  initial  model  disturbances.  It  is  shown 
that  for  a  rolling  model,  particular  care  must  be  exercised  in  reducing 
stability  data  from  free -flight  tests  in  a  wind  tunnel  wrhen  the  motion  is 
monitored  in  only  one  plane.  This  follows  as  the  orthogonal  components 
of  the  yawing  motion  of  a  rolling  model  can  differ  appreciably  and  can 
exhibit  pseudo- nonlinear  characteristics.  The  relationships  between  the 
frequency,  damping,  and  amplitude  of  the  motion  vectors  are  discussed 
in  relation  to  the  possible  variations  in  motion  patterns,  which  should  be 
an  aid  in  alleviating  some  apparent  misconceptions  that  have  arisen  from 
recent  free -flight  testing  in  wind  tunnels. 
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